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An example is provided on how a step can be digitized with , a high bandwidth utilizing heterodyning, 
it := .035 risetime of edge specified (ns) 
0 344 

f bw := f bw = 9.829 Bandwidth of critically damped second order system 

coO := 1 .554-2-rc-f bw calculate the center frequency for the system 



— = 1 5.274 center frequency (GHz) 

2-ji 



TD := 5 time delay for step edge (ns) 



2 

H(s) = — — e 



s + s + coO Is 

I Q ) 



laplace transform of the step specified 



the inverse Laplace transform provides the time-domain step waveform 
m if[t < TD,0,[-1 - <*Ht - TD)]-c [ - "° (t - TO)1 + l] 

To simulate the behavior of the analog components, it is modelled digitally with an extremely high sample rate 
FSj^ := 1 000 sample rate for simulating analog system (GHz) 

KH := 10000 kh := 0..KH - 1 

th. . := time of each point (ns) 

utilize a raised cosine window to minimize effects of the FFT 



wh. , := 
kh 



1 1 kh ^ 

= cod 2-7i 

2 2 | KH-1 ) 



wh = 1 1 thiS ° an be enab,ed to disab,e the windowing - essentially the same results are generated, 
w kh but the spectrums are not as pretty as with windowing. 

FIG. 4A 



calculate the windowed step 



.53181160838961202015 
t 10 := — + TD t lft = 5.006 



coO 



M0 : 



3 .889720 1 698674290579 

ton :== ~ — ■ ^ TD 

90 (DO Um^SMI 



t 90 - t|0 = 0.035 verify that risetime is correct 




Xh := CFFT(xh) Calculate the FFT 



NH:=^ nh:=0..NH fh •= J*.!^ 
2 nh NH 2 



step frequency content 




frequency (GHz) 



FIG. 4B 
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As we know, the scope does not have the bandwidth to digitize this signal. Therefore, we apply the method 
of this invention. First, we will utilize a system bandwidth of 5 GHz. then, we develop bandpass filters that 
select 5 GHz bands of the signal. Note that because the system is bandlimited, it is not actually necessary 
to utilize bandpass filters - only high pass filters need be utilized, but bandpass filters are used to simplify 
the discussion. Furthermore, the first band does not even need a filter - the scopes limited bandwidth will do 
this for us. (inside the scope, a digital low pass filter would be utilized to provide the hard bandwidth limiting) 

BW := 5 system bandwidth utilized for each band (GHz) 



make low pass and bandpass filters for each band 
nn := 1 .. NH - i 

Mfl nh := if Kh * BW > 1 - °) Mfh „h : = if ( BW < ^nh * 2 BW ' 1 .°) 
Mf Wnn Mfl NH-nn Mf Wnn : = Mfh NH- n „ 



Mfhh^ := if(2BW <fh nh < 3 BW, 1 ,u) 
Mft Wnn- Mf %H-„„ 
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Apply tnese niters to the input waverorm 
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frequency (GHz) 
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Xfl:=(Xh Mfl) 



Xfh:= (Xh Mfh) 



Xfhh := (Xh Mfhh) 



25 



25 



30 




30 
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Low Frequency Band 




High Frequency Band 




Very High Frequency Band 




15 

ftequehcy(GHz) 



xfl:= ICFFt(Xfi) 



xfh := ICTFTXXfh) xfhh := ICFFT(Xfhh) 



FIG. 4D 
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Here are the time domain waveforms at the output of the fitters 

Actual Input Waveform 




FIG. 4E 
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II is useful to add these three signals together and compare them to the input waveform. You will note the sum 
is not identical to the input because the system has limited the bandwidth at 15 GHz. The 15 GHz bandwidth 
limited signal is the best that we will be able to provide. 



Comparison 




mmmam Actual Input Step 

• • • Sum of Three Components 

It is also useful to compare the low frequency and actual input waveforms directly: 



Comparison 




4 5 6 

time (ns) 

"""" Actual Input Step 

• • • 5 GHz Bandwidth Limited Step 

The point of this last comparison is to demonstrate the problem that this invention is designed to solve. The limited 
bandwidth slows the edge of the step. This simulates the analog waveform that gets sampled by a digitizer with a 
front-end bandwidth of 5 GHz. Our goal is to digitize thejactual waveform with a much higher bandwidth. 

First, the high frequency and very high frequency bands are applied to the mixers 

F mixer0 := Bw *mixer0 :== md ( 2 n ) Tne frequency of the high frequency mixer is at the Cutoff frequency of the firs 

band. The frequency of the very high frequency mixer is twice that. 

F mixerl :== 2 BW *mixerl : = ™d(2*) 



FIG. 4F 
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apply the mixers 

xfhm kh := ^ft kh 2 cos(2- W .F mixel0 th kh + * mixtT0 ) *fi*"\h •= xfhK kh 2 co < 2 K F »nix C rl *jc h + *mixerl) 
look at the frequency content 

Xfhm := CFFT(xflim) Xfhhm := CFFT(xfhhm) 

Low pass filter the mixer outputs 



Xfhml := (XfhnvMfl) 



Xfhhm! := (Xflihm-Mfl) 



Note again that the typical manner of low pass filtering the mixer outputs would be to use the scope front-end. This 
filtering is being shown here as actual low pass filters applied. 



High Frequency Mixer Output 
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Very High Frequency Mixer Output 
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take the inverse FFT to generate the analog mixer output signals - the analog signals input to the 
channel digitizers. 



xfhml := ICFFT(Xfhml) xfhhml := ICFFT(Xfhhml) 

Low Frequency Waveform 




High Frequency Waveform 




I 

5 

time (ns) 



Very High Frequency Waveform 




FIG. 4H 



10/693,188 



it is interesting to see what the sum of these three waveforms are - there sums to not produce anything good 



Sum of Three Analog Waveforms 




At this point, the waveforms are digitized. The waveforms must be sampled at a rate sufficient to satisfy Nyquist 
For this example, this means that they must be sampled at at least 2 times BW. or 1 0 GS/s. After the waveforms 
have been digitized, they are immediately upsampled using SinX/x interpolation. This is possible because all 
digitized waveforms are bandlimited. It is useful to upsample the waveforms to a sample rate capable of meeting 
Nyquist for the system bandwidth - 1 have chosen 40 GS/s. The upsampling is trivial and for the purpose of this 
example, I simply use a 40 GS/s digitizer with the understanding that the exact same waveform would result 
from sampling the waveform at 10 GS/s and upsampling by a factor of 4. 

FS := 40 upsampled digitizer sample rate 
FS hi 

D : = — — D = 25 upsampling factor for analog waveform model 
rS 

K:=^ k:=0„K-l 
D 

sample the waveforms 



k FS 



\ : = X %D V =xfhml |c. D V =xh l 



k D 



W k := Wh k-D 



x hh := xfhhml k D 



generally, at this point, we would apply the sharp cutoff filter. If a sharp cutoff analog filter was not used we'd 
have to satisfy Nyquist such that any extra frequency content would not fold back Into the 5 GHz band. I've 
already applied a sharp cutoff filter to the analog signal, so this is not necessary. 

Also, at this point, some magnitude and phase compensation would probably be necessary to account for non- 
.deal channel frequency response characteristics. This example shows the signal digitized with ideal digitizers 
with ideal frequency response characteristics. 



FIG. 41 
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Next, the high and very nigh frequency waveforms are mixed up to there appropriate frequency location and 
digitally bandpass filtered. 

NOTE THAT THESE DIGITAL MIXERS KNOW THE PHASE OF THE ANALOG MIXERS - SOME MECHANISM 
MUST BE PROVIDED FOR DETERMINING THIS - EITHER THROUGH A PILOT TONE, OR LOCKING OF THE 
THE MIXER PHASE TO THE SAMPLE CLOCK. 

apply digital mixers 

*hm k - ^^^'^^^'^"^iiixcrO'S: "^~* ^roixcioj) *hhni k : = ^( 2 cos ( 2 jt ' F mixcrI ^ + *mixerl)) 

bandpass filter the mixer outputs 

N:= — n.-=0..N 
2 

f •= — — 
n N 2 

Xhm := CFFT(x hm ) Xhhm := CFFT(x hhm ) Xlm := CFFT(x,) 

Xfhm := ifff >BW,Xhm ,0\ Xfhhm := if If >2-BW,Xhhm ,0) 

nn:= I..N- 1 



Xfhm N + nn := Xfhm N-nn Xfhhm N + nn := Xfhhm N-nn 

X h := CFFT(x h ) X, := CFFT(x,) X hh := CFFT(x hh ) 



FIG. 4J 
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Low Frequency Waveform 
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High Frequency Digital Mixer Output 
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Very High Frequency Digital Mixer Output 
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Sum of AH Frequency Bands 



OO Sum 
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Actual Input 



By summing the output waveforms, we have acquired the waveform with a 15 GHz bandwidth utilizing three 
5 GHz bandwidth channels! 

Now lets see how the time domain waveforms compare 

The analog waveforms in the plots below are the analog outputs of the bandpass and low pass filters prior 
to the mixer 

xfhm := Re(ICFFT(Xfhm)) xfhhrn := Re(ICFFT(Xfhhm)) 
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Low Frequency Waveform 




™~ Analog Waveform 
O Digitized (and upsampled) Waveform 



High Frequency Waveform 




Analog Waveform 
O Digital Mixer and HP Filter Output 



Very H igh Frequency Waveform 
p— 
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Low Frequency Only 
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As you can see, the 15 GHz bandwidth limited step is recreated. 

Here are some risetime measurements: 
rt act :== riseTim^xh, FSjjj) 

rt Iow := riseTime(Re(xj),FS) 

rt high :== riseTiine(Re(xj) + Re(xfhm),FS) 

rtyhigh := riseTim^R^x,) + Rc(xfhm) + Re(xfhhm) f FS) 



bandwidth predicted bandwidth predicted 
using 0.35 multiplier using 0.45 multiplier 



35 =9.992 .45 



rt act 1 000 = 35.029 actual step risetime rt ~ y yyz j£L _ } 2 846 

act ^ 

rt Iow 1000 = 93.196 5 GHz bandwidth -^L = 3.755 A5 

risetime rt Iow — = 4.829 

"low 

= 8.183 



1 aaa caa 10 GHz bandwidth .35 
^•1000 = 54.99 risetime — = 6.365 ^45 

rt high 



rt high 
.45 

risetime rtyhjgh rtyhigh 



^high 1000 = 43 73 15 GHz bandwidth — = 8.004 45 = 10 .29 



rtv high 10 = 0.437 multiplier determined by 10 GHz bandwidth (noting that the signal 
itself only required 1 0 Ghjz of bandwidth) 



FIG. 4P 
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An example is provided on how a step can be digitized with a high bandwidth utilizing heterodyning, 
rt := .045 risetime of edge specified (ns) 

0 144 

f bw := r f bw = 7.644 Bandwidth of critically damped second order system 

toO := 1 .554-2 rc f bw calculate the center frequency for the system 

— = 1 1 .879 center frequency (GHz) 

2-71 

TD := 5 time delay for step edge (ns) 

2 

~ _ , v toO -*TD 

H(s) = — : — e 

laplace transform of the step specified 



( 2 caO 2\ 
s + s+ coO \s 



the inverse Laplace transform provides the time-domain step waveform 
f(t) := ifLt < TD,0,L-1 ~ co0 (t - TD)h l - <S30<i - TD)] + 1 J 

To simulate the behavior of the analog components, it is modelled digitally with an extremely high sample rate 
FS hi := 1000 sample rate for simulating analog system (GHz) 

KH:= 10000 kh:= 0..KH- 1 

th kh = ~~~ tJme °* each point ^ ns ^ 
FS hi 

utilize a raised cosine window to minimize effects of the FFT 



wh. . 
kh 



1 1 (n kh ^ 

:= cos 2 k 

2 2 ^ KH - 1 ) 



■ this can be enabled to disable the windowing - essentially the same results are generated, 
wh kh : ~ 1 but the spectrums are not as pretty as with windowing. 



FIG. 5A 
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Lh kh : " ^^i^kh calculate the windowed step 



.53181160838961202015 



10 



+ TD t 1Q = 5.007 



3.8897201698674290579 __ 
t 90 := + TD 



coO 



t 90 = 5.052 

t 90 - t 10 = 0.045 verify that risetime is correct 
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Xh := CFFT(xh) Calculate the FFT 
KH 



NH 



nh FS hi 



nh:= 0..NH fh . := 

nh ^ 2 



step frequency content 
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As we know, the scope does not have the bandwidth to digitize this signal. Therefore, we apply the method 
of this invention. First, we will utilize a system bandwidth of 5 GHz. then, we develop bandpass filters that 
select 5 GHz bands of the signal. Note that because the system is bandlimited, it is not actually necessary 
to utilize bandpass filters - only high pass filters need be utilized, but bandpass filters are used to simplify 
the discussion. Furthermore, the first band does not even need a filter - the scopes limited bandwidth will do 
this for us. (inside the scope, a digital low pass filter would be utilized to provide the hard bandwidth limiting) 



BW := 5 system bandwidth utilized for each band (GHz) 

make low pass and bandpass filters for each band 

nn:= 1..NH- 1 

Mfl nh := if Kh * BW ' 1 - °) Mfh nh : = if ( BW <fh *h* 2 BW > 1 - °) 
Mf Wnn := = Mf W-nn Mfh NH+nn Mfh NH-nn 



Low Frequency Band 



I 



10 



15 

frequency (GHz) 



20 



25 



30 



High Frequency Band 



-3 
B 



10 



15 

frequency (GHz) 



20 



25 



30 



Apply these filters to the input waveform 



Xfl := (Xh-Mfl) 



Xfh := (XhMfh) 
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xfl := ICFFT(Xfl) 



xfh := ICFFT(Xfh) 



FIG. 5D 
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Here are the time domain waveforms at the output of the filters 

Actual Input Waveform 
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It is useful to add these three signals together and compare them to the Input waveform. You will note the sum 
is not identical to the input because the system has limited the bandwidth at 15 GHz. The 15 GHz bandwidth 
limited signal is the best that we will be able to provide. 



Comparison 




Actual Input Step 
• • • Sum of Three Components 

It is also useful to compare the low frequency and actual input waveforms directly: 



Comparison 



r. 

is . 




5 

time (ns) 



Actual Input Step 

5 GHz Bandwidth Limited Step 



The point of this last comparison is to demonstrate the problem that this invention is designed to solve. The limited 
bandwidth slows the edge of the step. This simulates the analog waveform that gets sampled by a digitizer with a 
front-end bandwidth of 5 GHz. Our goal Is to digitize the actual waveform with a much higher bandwidth. 

First, the high frequency and very high frequency bands are applied to the mixers 

F mixerO := 2 Bw °mixer0 := rnd(2-rc) The frequency of the high frequency mixer is at the Cutoff frequency of the first 

band. 



FIG. 5F 
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apply the mixers 

xfhm kh := ^kh^ ^^^^mixeiO ^ + ^mixeio) 
look at the frequency content 
Xfhm := CFFT(xfhm) 
Low pass filter the mixer outputs 



Xfhml := (XfhnvMfl) 

Note again that the typical manner of low pass filtering the mixer outputs would be to use the scope front-end. This 
filtering is being shown here as actual low pass filters applied. 
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take the inverse FFT to generate the analog mixer output signals - the analog signals input to the 
channel digitizers. 





FIG. 5H 
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it is interesting to see what the sum of these three waveforms are - there sums to not produce anything good 



Sum of Three Analog Waveforms 




At this point, the waveforms are digitized. The waveforms must be sampled at a rate sufficient to satisfy Nyquist. 
For this example, this means that they must be sampled at at least 2 times BW, or 10 GS/s. After the waveforms 
have been digitized, they are immediately upsampled using SinX/x interpolation. This is possible because all 
digitized waveforms are bandlimited. It is useful to upsampfe the waveforms to a sample rate capable of meeting 
Nyquist for the system bandwidth - 1 have chosen 40 GS/s. The upsampling is trivial and for the purpose of this 
example, I simply use a 40 GS/s digitizer with the understanding that the exact same waveform would result 
from sampling the waveform at 10 GS/s and upsampling by a factor of 4. 



FS := 40 upsampled digitizer sample rate 
FS hi 



D:= 



FS 



D = 25 upsampling factor for analog waveform model 



K:=™ k: =o..K-l 
D 



sample the waveforms 



* k ~ FS 



\ := ^k-D X \ := Xmm ^k D \ := ^kl 



w k := wh k . D 



generally, at this point, we would apply the sharp cutoff filter. If a sharp cutoff analog filter was not used, we'd 
have to satisfy Nyquist such that any extra frequency content would not fold back into the 5 GHz band. I've 
already applied a sharp cutoff filter to the analog signal, so this is not necessary. 



Also, at this point, some magnitude and phase compensation would probably be necessary to account for non- 
ideal channel frequency response characteristics. This example shows the signal digitized with ideal digitizers 
with ideal frequency response characteristics. 
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Next, the high and very high frequency waveforms are mixed up to there appropriate frequency location and 
digitally bandpass filtered. 

NOTE THAT THESE DIGITAL MIXERS KNOW THE PHASE OF THE ANALOG MIXERS - SOME MECHANISM 
MUST BE PROVIDED FOR DETERMINING THIS - EITHER THROUGH A PILOT TONE, OR LOCKING OF THE 
THE MIXER PHASE TO THE SAMPLE CLOCK. 

apply digital mixers 

x hm k : = Xt 2 ' 008 ! 2 '"' 17 ^ 10 ^ + ^mixero)) 
bandpass filter the mixer outputs 

N:= — n:=0..N 



2 




FS 



2 




Xlm := CFFT(xj) 



nn:= 1..N - 



Xfhm N+nn := Xfhm N-i 
X h := CFFT(x h ) 



•nn 



X! := CFFT(x 1 ) 
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Low Frequency Waveform 
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Sum of All Frequency Bands 
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By summing the output waveforms, we have acquired the waveform with a 15 GHz bandwidth utilizing three 
5 GHz bandwidth channels! 



Now lets see how the time domain waveforms compare 

The analog waveforms in the plots below are the analog outputs of the bandpass and low pass filters prior 
to the mixer 



xfhm := Re(ICFFT(Xfhm)) 
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As you can see, the 15 GHz bandwidth limited step is recreated. 

Here are some risetime measurements: 
rt act ::= riseTime(xh,FS hi ) 

rt, ow := riscTime(Re(x,),FS) 

rt high := nseTime(Re(x|) + Re(xfhm),FS) 



bandwidth predicted bandwidth predicted 
using 0.35 multiplier using 0.45 multiplier 



35 

J3 = 7.772 .45 



rt act - 1000 = 45.036 actual step risetime rtact - — = 9.992 

30 rt act 



rt low 1000 = 95.195 5 GHz bandwidth = 3.677 .45 



risetime rt l< 



low rtl 



4.727 



low 



• /w* rt ****** 10 GHz bandwidth .35 Ac 

rt high - 1000 = 59.706 — — = 5.862 _A5_ = 

^ rt hi gh 

rthigh* 10 = 0 597 multiplier determined by 10 GHz bandwidth (noting that the signal 
itself only required 1 0 GHz of bandwidth) 
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